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^ ■ Abstract 

' . IGRJ1736 1-4441 is a hard transient recently observed by the INTEGRAL 

^ \ satelhte. The source, close to the center of gravity of the globular cluster 

tg \ NGC 6388, quickly became the target of follow-up observations conducted 

by the Chandra, Swift/XRT and RXTE observatories. Here, we concentrate 



(N 



in particular on a set of observations conducted by the XMM- Newton satel- 
> . lite during two slews, in order to get the spectral information of the source 

/O \ and search for spectral variations. The spectral parameters determined by 

Q\ \ the recent XMM-Newton slew observations were compared to the previously 

^ ' known results. The maximum unabsorbed X-ray flux in the 0.5-10 keV band 

^ ■ as detected by the XMM-Newton slew observations is ~ 4.5 x 10^"*^^ erg cm~^ 

^ . s~^, i.e. consistent with that observed by the Swift/XRT satellite 15 days 

earlier. The spectrum seems to be marginally consistent (F ~ 0.93 — 1.63) 
with that derived from the previous high energy observations. 

rN ■ Keywords: (Galaxy:) globular clusters: individual (NGC 6388), X-rays: 

c3 . general 



1. Introduction 



IGRJ17361-4441 is a hard X-ray tra nsient source first obs erved by lGibaud et al. 



( 1201 if l with the IBIS /ISGRI telescope flUbertini et al.l .l2003l) onboard the IN- 



TEGRAL satellite ( 1 Winkler et al.l . |2003[ ) and was quickly recognized t o be 



hosted in the galactic globular cluster NGC 6388 (JFerrigno et al.l . 1201 ll ). 
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The location of the transient (close to the globular cluster gravitational 
center, but see later) is of great importance since NGC 63 88, among all the 



globu lar clusters in our Galaxy, is one of the best candidates (iBaumgardt et al. 



20051 ) to host an intermediate mass black hole (he reafter IMBH) . In pa rticu- 
lar, by using high resolution optical observations, iLanzoni et al.l (120071 ) esti- 
mated the mass of the IMBH to be ~ 5700 Mq. It would be natural for such 
an IMBH to emit significant radiation in the X-ray band due to the likely 
accretion of matter from its su rroundings. In the context of t he earliest ob- 
servat ions of globular clusters, iBahcall fc Ostrikerl (119751 ) and ISilk fc Arons 
(119761 ) were the first to suggest that the X-ray emission detected towards 



these clusters was due to IMBHs (in the mass range 20 Mq - IO^Mq) accret- 
in g from the intracluste r medium. This issue was considered more recently 
by lGrindlay et al.l (120011 ) who provided the census of the compact object and 
binary population in the globular cluster 47 Tuc and obtained an upper limit 
to the central IMBH of a few hundred solar masses. 



Initial XMM- Newton an d Chandra observations in this direction (iNucita et al, 

20081 and lCsehet"aDl2010h showed that the core of NGC 6388 hosts several 
X-ray sources. B ased on the co r relati on between the X-ra y and radio flux 
from black holes (JMerloni et al.l (120031 )): iMaccarond (120041 ) was the first to 
point out that the search for radio emission from faint black holes is useful to 
test the IMBH hypothe sis in globular cluste rs and dwarf spheroidal galaxies 
(JMaccarone et al.l |2005|) . ICseh et al.l (120101 ) observed the central region of 
NGC 6388 in the radio band using the Australia Telescope Compact Array 
(ATCA) to search for radio signatures of the IMBH. The radio observation 
resulted in an upper limit of the IMBH mass of ~ 1500 Mq. 

The discovery of a transient source close to the NGC 6388 gravitational 
center could be related to the turning on of the putative globular cluster 
IMBH. However, as will become clear in the subsequent sections, the nature 
and spectral properties of the transient IGRJ17361-4441 are difficult to rec- 
oncile with the IMBH picture and rather favour an interpretation as a high 
mass X-ray binary (HMXB) or a low mass X-ray binary (LMXB). Several 
observational campaigns (in the X-rays as well as in the radio band) were 
organized in order to pinpoint IGRJ17361-4441 and draw firm conclusions 
on the NGC 6388 IMBH paradigm. 

In this paper we briefiy discuss (see se ction |2l) the past X -ray o bservations 



of th e NGC 6388 globular cluster (see iNucita et al.l I2008L and ICseh et al 
2010h a nd the discovery of the hard transient IGRJ17361-4441 by INTE- 
GRAL (JGibaud et al.ll201ll ) as well as the follow-up observations conducted 



by Ch an dra (Poolev et al.ll2011l) . Swift /XRT and RXTE observatories ( Ferrigno et al 
201ll and lBozzo et al.ll201ll ). Then we concentrate (see section[3]) on the anal- 
syis of two XMM-Newton slew observations of NGC 6388 conducted 15 days 
after the INTEGRAL discovery of the source. The two slew observations 
had ~ 7.6 seconds and ~ 7.7 seconds on source exposure time. Finally our 
conlusions are presented in section |H 



2. Previous observations of NGC 6388 



2.1. XMM-Newton and Chandra observations of NGC 6388 

By studying a combination of high resolution (HST ACS-HRC, ACS- 
WFC, and WFPC2 ) and wide fiel d (ESO -WFI) observations of the globular 
cluster NGC 6388, iLanzoni et al.l (120071 ) claimed the existence of a central 
IMBH. Such a compact object of mass ~ 5700 M© should reside in the 
globular cluster center of gravit}u localized at the coordinates ( J2000) RA = 
17^ 3gm i7_23«, Dec = -44° 44' 7.1". An uncertainty of 0.3" is associated 
with both coordinates. 



Nucita et al.l ( l2008l ) suggested that this IMBH should emit radiation in 



the X-ray band due to accretion from the surrounding matter. A 48 ks XMM- 
Newton observation was made on 21 March 2003. It resulted in a spectrum 
which was well fit by an absorbed power-law model. The resulting best fit 
parameters were N^ = (2.7 ± 0.3) x 10^^ cm~^ for the hydrogen column 
density and F = 2.4 ±0.1 for the power law index. The unabsorbed flux in 
the 0.5 — 7 keV band was -Fo.5_7 = (4.0 ±0.2) x 10~^^ erg cm"^ s~^ which, for a 
distance of 13.2 kpc corresponds to a luminosity of Lo.5-7 — (7.2 ±0.4) x 10^^ 
erg s~^. Note that the hydrogen column density is consistent with t he average 



one found in the direction of the target ( Dickey fc Lockmanlll990l ) 



^Note however that, as first pointed out bv iBahcall fc Wolj ( 19761 ). a black hole in a 
stellar cluster will experience a Brownian motion due to gravitational interactions with 
the surrounding objects. Thus, the black hole is not necessarily at th e dynamical center o f 
the host cluster, but may move with mean square velocity given by ([Merritt et al.ll2007() 



1 3 



(1) 



where M, m and a represent the black hole mass, the perturber average mass and the 
stellar velocity dispersion wit hin ~ 0.6ri, respec tively. Here r^ is the influence radius of 
the black hole (for details see iMerritt et al.ll2007l and references therein). 



The Chandra satelhte, with a much better angular resolution than that 
of XMM- Newton, observed towards NGC 6388 for ~ 45 ks on 21 April 2005 



(id 5505). iNucita et al.l ( 20081) identified 16 discrete sources within the half 
mass radius (~ 40", see iLanzoni et al.l 120071 ) of the cluster. The 3 sources 



close to the gravitational center were not spatially resolved by the authors, so 
that they were considered virtually as a single source (labeled as #14*). The 
unabsorbed flux in the 0.5-7 keV band of the #14* is -Fo.5_7 ~ 1.7 x 10"^^ 
erg cm~^ s~^, corresponding to a luminosity of -Lo.5-7 ~ 3 x 10^'^ erg s~^. 

A more de t ailed analysis on the same Chandra data set was conducted 
by lCseh et al.l ( 20101 ). After removing the pixel randomization, these authors 
were able to spatially resolve the source #14* into three separate sources 
labeled as #12, #7 and #3. In particular, the source #12, which is consistent 
with the position of the center of gravity of NGC 6388, is characterized by 



an unabsorbed flux of Fq.; 
int rinsic lumiriqsiy of L, 



-14 



0.5-7 



'. 4.0 X 10 "■" erg cm" 
8.3 X 10^2 grg s-^ 



s corresponding to an 



Cseh et al.l (120 lOf ) searched for a radio counterpart of the putative IMBH 
in NGC 6388 using the ATCA facility. Unfortunately, this search only re- 
sulted in an upper limit to the radio flux at 5 GHz of ~ 81 /iJy/beam. 
Therefore, it was only possible to determine an upper limit to the IMBH 
radio luminosity of Lr < 8.4 x 10^^ erg s~^. 



Based on the fundaniental plane of black hole accretion (iMerloni et al. 



2OO3I and iKoerding et al.l |2006| ) and using the observed X-ray and radio lu- 



minosities, it was then possible to p ut a 3a uppe r limit of ~ 1500 Mq on 



the mass of the IMBH in NGC 6388 (ICseh et al.ll20ia) . The estimated mass 
value has to be treated with caution for two reasons: i) the identification of 
the X-ray counterpart of such a black hole is not trivial since several sources 
are close to the NGC 6388 center of gravity. If none of them are associ- 
ated with the IMBH, then one can not use the fundamental plane relation to 
ge t an estimate of the ma ss; ii) the fundamenta l plane relation (as derived 
by iMerloni et al.l |2003| and iKoerding et al.l |2006|) is not tested for black hole 
masses in the range of inte r est fo r IM BHs, i.e. 10 ^ M^^ - 10^ M^^. Note how- 



ia. 



Hi 



ever that iMaccarone et al.l (120051 ) and iMaccarone fc ServillatI (120081 ) showed 
that the non-detection of a radio source, in combination with the estimate 
of the globular cluster ISM densitjo and the expected value of the accre- 



•^The amount of gas contained in globular clusters is an issue of debate. The intracluster 
medium density can be estimated by using the dispersion measures of the pulsars observed 



tion rate, can be used to get information (at least as an order of magnitude 
estimate) of the IMBH mass. 

2.2. INTEGRAL discovery of IGRJ17361-4441 and subsequent X-ray follow- 
up observations 



On 11 August 2011, iGibaud et al.l ( 120111 ) reported the discovery of a 
new hard X-ray tran sient (IGRJ17361-4441) by the IB IS/ISGRI telescope 
flUbertini et al.l . l2003h onboard the INTEGRAL satellite dWinkler et al.l . l2003h . 
The spectrum of the source, associated with the globular cluster NGC 6388, 
was described by a power law with photon index T 
ized by a flux in the 20-100 keV of -F20-100 — 

Since this newly discovered transient is possibly associated with the IMBH 
in NGC 6388, IGRJ17361-4441 became the target of several X-ray follow-up 
observations aimed at obtaining accurate position and flux measurements to 
tes t if this associat i on is correct. 

Wijnands et al.l ( 1201 ll ) (but see also lFerrigno et al.ll201ll ) reported a Swift /XRT 



2.61q;7 ^^d character- 
9.7 X 10"^^ erg cm~^ s" 



observation (1.9 Ks on 16 August 2011) in which the astrometrically corrected 
position of IGRJ17361-4441 was RA = 17^ 36°^ 17.5^ Dec = -44° 44' 7.1". 



A more detailed analysis on the XRT data conducted by lBozzo et al.l (120111 ) 
determined the new transient position to be RA = 17*^ 36™ 17.27**, Dec = 
—44° 44' 7.0" with an associated error (on both coordinates) of 1.9". Thus, 
the distancqj of the transient from the center of gravity of NGC 6388 is 0.4" ± 
1.4". Hence, the source position is consistent (according to the Swift /XRT 
data) with the center of gravity of the cluster and (possibly) associated with 
the IMBH. Note also that a 2.5 Ks Chandra observation was made on 29 



August 2011 (jWijnands et al.l 1201 ll ) in order to improve the accuracy of 
the location of the transient. IGRJ17361-4441 is located at the coordinates 
RA = 17'^ 36™ 17.418% Dec = -44° 44' 5.98" (the nominal Chandra po- 
sitional accuracy is 0.6" on both coordinates). In this case, the estimated 
distance of the transient to the cluster center of gravity is 2.3" ± 0.5". 

Based on the Swift/XRT astrometry only, one could conclude that the 
position of IGRJ17361-4441 is formally consistent with the center of gravity 



within the cluster (see e.g. iFreire et aLll200lL 120031) or inferred by the empirical knowledge 
about the stellar mass loss (jPfahl fc Rappaportll200l[ ). 

■^ We note that the distances between the sources and the center of gravity of the globular 
cluser NGC 6388 are calculated by using the well know Haversine formula. Distance 
uncertainties are calculated by corrctly propagating the errors on both a and S coordinates. 



of NGC 6388 and possibly related to the putative IMBH in the globular 



cluste r. An updated radio observation conducted at ATCA by iBozzo et al. 



(120 111 ) put a more stringent upper limit to the radio lum inosity of Lr < 



5 X 10 erg s so that, following the same procedure as in lCseh et al.ll2010 



and the 2005 Chandr a X-ray flux estini ate, the new IMBH upper limit turns 
out to be ~ 600 Mg (IBozzo et al.ll201lh . 

However, a caveat on this conclusion is necessary. The new Chandra re- 
fined source coordinates (even if formally consistent with the source position 
determined by Swift/XRT) indicate that the transient could be a new X- 
ray source (see later) not associated with the IMBH. In this case, and for 
the reasons explained above, one should not use the black hole fundamental 
plane relation in order to estimate the IMBH mass. If one believes that the 
transient is associated with the NGC 6388 center of gravity, then it should 
also be noted that at least three sources (those labeled as #12, #7 and 
#3 in lCseh et al.ll2010l ) are within the error box of Swift/XRT. In particular, 
sources #12 and #7 have fluxes ^ 4.0 x 10~^^ erg cm~^ s~^ and ^ 6.9 x 10^^"^ 
erg cm~^ s~^, respectively. If source 7 is associated with the IMBH, the X- 
ray and radio observations together with the fundamental plane relation give 
an upper limit of ~ 1200 Mq. 



The Swift/XRT spectrum of the transient source was fitted (IBozzo et al. 



20111 ) with an absorbed power law with photon index F ~ 0.5 — 0.9 and 



hydrogen column density Nh — (0.5 — 0.9) x 10^^ cm ^, i.e. consistent with 
that derived from the previous XMM-Newton data. 

The flux in the 1.0-10 keV band is Fi_io = (4.5 - 4.8) x 10" 



11 



erg cm 



'^ a factor 100 more l uminous than the source #12 possibly associated with 

2010) ■ When the broad band spectrum 



the NGC 6388 IMBH (iCseh et al 



(obtained by using Swift/XR T and INTEGRAL /ISGRI) was analysed (and 
fit with a broken power law), IBozzo et al.l ( 201ll ) obtained a flux of -Fi_io = 
(4.6^q;5) X 10~^^ erg cm"^ s~^ in the 1.0 — 10 kev band and a flux of -F20-100 = 
(7.8+3I) X 10-^1 erg cm"2 g-^ in the 20 - 100 kev band. 

These results are consistent with those obtained by using the XRTE/PCA 
follow-up observation made on 17 August 2011. In particular, i t was found 
that Fs- 15 = (6.7+11:4) X 10"^^ erg cm'^ s"^ in the 3-15 keV band (IBozzo et all 
201lh . 



3. The XMM-Newton slew observations 



The large collecting area of the nested mirrors together with the high 
quantum efficiency of the EPIC/PN camera make the XMM- Newton satel- 
lite t he most sensitive X-ray observatory available at present (jJansen et al. 
200ll ). The XMM-Newton satellite was recognized to be a good instrument 
to collect d ata during slewing m anoeuvres and is performing an X-ray survey 
of the sky (jSaxton et al.ll2008l ). Note that the XMM-Newton slew observa- 



tions are moderately deep with a detection limit of 1.2 x 10 ^^ erg cm ^ s ^ 
in the 0.2 - 12.0 keV band. 

Due to the scheduled observation program, the XMM-Newton satellite 
twice observed the region of the sky around NGC 6388. The observations 
were taken on September 1'^* 2011 at 13:10:34 (herafter SI observation) 
and 19:00:17 (UT) (herafter S2 observation), i.e. 15 days after the first 
Swift/XRT follow-up observation of IGRJ17361-4441. The transient source 
was then observed serendipitousljQ for ~ 7.6 s and ~ 7.7 s in the two slew 
observations with the EPIC/PN instrument (see Fig. [T]). 

We decided to analyze the data sets separately in order to study a possible 
spectral variation on time-scales of hours. 

The source spectrum has been extracted from a circle of radius 60" about 
the source with the background being extracted from an annulus of inner ra- 
dius 90" and outer radius 120" about the source. The detector matrices are 
calculated taking into account th e transit of t he sou rce across the detector 
and using the method described in lRead et al.l (120081 ). Hence, the source and 
background spectra (as well as the response matrices) were imported in the 
XSPEC package (version 12.4.0) for the spectral analysis and fitting proce- 
dure. The adopted model is an absorbed power law {wabs*power) with the 
hydrogen column density fixed to the average value found by ROSAT in the 
direction of the target, i.e. 2.5 x 10^^ cm~^. Note that this value is consistent 



with that derived by iNucita et al.l (120081 ) when analyzing the 2005 XMM- 
Newton observation of NG C 6388 and also similar to t he column density 
found by lBozzo et al.l ( 20111 ) (see also lFerrigno et al.ll201ll ) while studying the 
Swift/XRT follow-up observation of IGRJ17361-4441. The adopted model 



'*The exposure times were estimated by calculating the distance travelled by the source 
on the detector with the typical XMM-Newton slew speed of 90 deg h~^. T he resulting 
exposure times were also corrected for chip gaps and bad pixels (for details see lRead et al 
[2008ii . 
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Figure 1: Contours (increasing by factors of 2) of lightly-smoothed 0.2-12 keV XMM- 
Newton slew emission (EPIC-pn camera, the two 01/09/11 observations combined), su- 
perimposed on a SAO-DSS image of NGC6388. 
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Figure 2: The XMM-Newton spectrum of IGRJ17361-4441 (data points) collected during 
the first slew observation on September F** 2011. The solid line represents the best fit 
model (see text for details). 
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Figure 3: The same as in figure 1, but for the second XMM-Newton slew observation on 
September I'** 2011. 



has two free parameters; the photon index F and the power law normahzation 
N. 

The fitting procedure to the SI spectrum resuhed in the best fit param- 
eters {x^/u = 1.4 for 11 d.o.f.) r = 1.16 ± 0.20 and A^ = (1.7 ± 0.7) x 10~l 
The absorbed fluxes in the 0.5 - 2.0 keV, 2.0 - 10 keV, and 0.5 - 10.0 keV 
bands are Fo.5_2 = (5.4+2-8) x lO'^^ g^g cm^s g-i^ p^_^^ ^ (3.6+};^) x lO^^^ 
erg cm~2 s~\ and -Fo.s-io = (4.1+};9) x 10"^-*^ erg cm"^ s~^, respectively. 

Fitting the S2 spectrum with the same model gives the best fit parameters 
(x^/u = 0.4 for 10 d.o.f.) F = 1.28 ± 0.35 and A^ = (1.2 ± 0.6) x lO'^. The 
absorbed fluxes in the same bands as above are -Fo.5_2 = (3.4 



-12 



1.9+1:4) X 10~^2 erg cm-2 s-\ and Ff 



0.5-10 



+li) X 10"^- erg 

+1.2\ 



(2.1++: 



X 



cm ' s % i^2-io 

10~^^ erg cm^2 s^^, respectively. 

The maximum unabsorbed flux of the source in the 0.5-10 keV band 
is ~ 4.5 X 10^^^ erg cm"^ s~^ corresponding to an intrinsic luminosity of 
~ 9.3 X 10^^ erg s"^ 

In Figure HI the unabsorbed X-ray fluxes of NGC 6388 in the 0.5 — 10 
keV band (from 2003 to 2011) are shown. In the insert, we give the data 
points corresponding to the IGRJ17361-4441 flare observed and monitored by 
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several instruments (Swift/XRT, RXTE/PCA and XMM-Newton) in 2011 
only. Note that for the XMM-Newton observation in 2003, the Chandra 
observation in 2005, and the Swift and RXTE in 2011 the 0.5 - 10 keV 
band fluxes were obtained by extrapolating fto this energy ba nd) the best 
fit models available in the literature (see e.g. iBozzo et al.ll201ll ). 

What can be easily observed from the results reported here, is that the 
X-ray flux in the 0.5-10 keV band as detected by the XMM-Newton slew 
observations, i.e. ~ 4.5 x 10~^^ erg cm~^ s~^, is consistent with that observed 
by Swift/XRT 15 days earlier. Note also that the spectrum power-law seems 
to be marginally consistent (F ~ 0.93 — 1.63) with th at derived from th e 
previous high energy observations (F c^ 0.5 — 0.9, see IBozzo et al.l 120111 ). 
a point that could help in a classification of the transient nature (see the 
subsequent discussion) . 
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Figure 4: The data points from the left to the right correspond to the NGC 6388 flux in the 
0.5 — 10 keV band from 2003 to 2011. In the insert, we give the data points corresponding 
to the flare observed and monitored in 2011 only (see text for details). 



4. Results and discussion 

IGRJ17361-4441 is a hard transient recently observed by the INTEGRAL 
satellite. X-ray follow-up observations have shown that the source is within 
the globular cluster NGC 6388. Based only on the astrometry of the Swift/XRT 
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satellite, the transient position is consistent with the center of gravity of the 
globular cluster and this opens the possibility that IGRJ17361-4441 is asso- 
ciated with an IMBH which is turning on. 

However, if one believes that the transient is associated with the IMBH in 
NGC 6388, then it should be n oted that at leas t three X-ray sources (those 
labeled as #12, #7 and #3 in lCseh et al.ll2010[ ) are within the error box of 
Swift/XRT. In particular, the sources #12 and #7 have fluxes which differ 
at least by a factor 2 between them. If source #7 is associated with the 



IMBH, then the observed Chandra X-ray flux (see iNucita et al.l l2008l and 
Cseh et al.l 120101 ) and the updated radio observation of IBozzo et al.l (120111 ) 
together with the fundamental plane relation give an upper limit of ~ 1200 

Mq. 

In the IMBH hypothesis, the intrinsic luminosity of the source as deter- 
mined by using the XMM-Newton slew data (i.e. ~ 9.3 x 10^^ erg s~^) 
should be compared with that de rived by using the 2005 Chandra data (and 
in particular for the source #12 in lCseh et al.ll2010[ ) when the putative IMBH 
was in quiescent state. In this case, one finds that the source luminosity in- 
creased by at least a factor ~ 1000. Moreover, the spectrum seems to follow 
a power law with photon index F ~ 0.96 — 1.63. 

Nevertheless, the refined source position given by the Chandra satellite 
(even if still in agreement with the Swift/XRT result) argues against the 
IMBH hypothesis in favor of a newly discovered source. In this case, the 
XMM-Newton intrinsic source luminosity sh ould be compared w' ith the up- 
per limit for the quiescent state of the source. iPooley et al.l (l201ll ). based on 
the non-detection of the source in the 2005 Chandra observation, estimated 
this limit to be ~ lO^^erg s^^. Thus, in this case the transient source has 
increased its luminosity by a factor close to 10^. 

Two more possibilities for the nature of the transient source are that 
it is either a HMXB or a LMXB. The first possibility is actually unlikely 
since these systems involve c ompanion stars with mass larger than ~ 10 M© 
( Lewin &: van der Klisll2006l ). i.e. 0/B stars which are not expected to exist 
in globular clusters. Note also that NGC 6388 was extensively observed by 



the HST instruments (see e.g. ICatelan et al.l 120061 ) and the collected data 



did not show the presence of any O or B star in the globular cluster. 

Hence, among the IMBH alternatives, the LMXB option is the most 
favorable. This is supported by the X-ray luminosity (^ 9.3 x 10'^^ erg s~^) 
and by the soft spectrum (F ^ 0.93 — 1.63) observed in the XMM-Newton 
slew observation which seems to be consistent with the typical characteristics 
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of the LMXB class of objects. A long X-ray observation (sufficient to allow a 
detailed timing and spectral analysis) may help in understanding the physics 
underlying this transient source. 
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